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Abstract 

A SUSY SO(IO) Xj44 gut model is constructed for fermion masses and mixing by 
introducing a minimal set of low dimensional Higgs representations needed to break 
the guage symmetry down to SU{3)c x U{l)em- The hierarchy of fermion masses can be 
^ understood in the framework of A4 symmetry. From the A4-invariant superpotential, 

the "double lopsided" mass matrices for the charged leptons and the down quarks are 
obtained. It is shown that this structure leads to bi-large neutrino mixings simultane- 
QJ ' ously with small CKM mixing angles. An excellent fit to the masses and mixings of 

the quarks and leptons as well as to CP violation parameter is obtained. Moreover, 
the model predicts the neutrino mixing angle sin^is ~ 0.15. 

(N 
> 

\o ; 1 Introduction 

q{ . There are many unification models that try to combine the strong and electroweak inter- 
O , 

0\ . actions into a simple group. The simplest model is based on SU(5) gauge symmetry [1]. 
O 

The minimal SU(5) model predicts the good mass relation for the third generation (i.e. 
^ ' = m°) at GUT scale. However, it has a bad prediction for mass relations for the first 

H : 

. . . and second generations (i.e. m° = m°, m° = m^) at GUT scale. In addition, it does not 
naturally accommodate the right handed neutrino. On the other hand, the SO (10) model 
accommodates all chiral fermions of one standard model generation plus a right handed 
neutrino within a 16-dimensional irreducible representation (irrep.). Also, minimal SO (10) 
with only IOh leads to the up quark mass matrix being proportional to the down quark 
mass matrix, so it is considered a good zeroth order approximation for CKM mixings. Some 
models based on SO (10) symmetry without including any family symmetry were proposed 
to explain most of the features of quarks and leptons However, one is not really 
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fully satisfied with only producing the fermion masses and mixing angles without explaining 
why we have three generations and without understanding the relation among generations, 
such as hierarchy and mixing angles. For example, the flavor symmetry can be 

employed to explain why the observed neutrino mixing matrix is in very good agreement 
with the so called tri-bi-maximal(TBM) mixing structure [6]. Thus, it may be important to 
consider the underlying family symmetry. One of the best candidates for flavor symmetry is 
the non-Abelian discrete symmetry A^, for the following reasons. Firstly, it is the smallest 
group that has a 3-dimensional irrep. Secondly, SUSY-SO(IO) x symmetry solves FCNC 
problem since the scalar fermions, which belong to the 16-irrep of SO(IO) and transform as 
a triplet under A4, have degenerate masses. Finally, it was shown that the TBM mixing 
structure can be obtained by imposing A4 symmetry [5]. 

Few models based on the SO(IO) x A4 group have been studied [7][H]- In these models, 
large Higgs representations have been employed. For example, in ref. [7], the authors em- 
ployed a (126//, 3) representation, where the first (second) entry indicates the transformation 
under SO(IO) (^44), in order to produce the fermion masses and mixing angles for both nor- 
mal and inverted neutrino mass spectra. Besides employing the large Higgs representation 
126//, the models in ref. [8] contain more than one adjoint 45// representation. It has been 
shown that only one adjoint Higgs field is required to break SO (10) while preserving the 
gauge coupling unification [9j. Also, using quite large Higgs representation like 126// leads 
to the unified gauge coupling being non-perturbative before the Planck scale, which might 
be hard to obtain from superstring theory [10]. Therefore, the purpose of this paper is to 
construct SO(IO) x ^4 model in which SO(IO) is broken to the standard model (SM) group 
in the minimal breaking scheme. This means using only a spinor-antispinor (16//, 16//) to 
break the rank of SO(IO) from five to four, and the right handed neutrino gets a heavy 
mass from the antispinor Higgs field (16//). Then one adjoint representation 45// is used 
to break the group all the way to the SM group. Recently, a numerical analysis for quark 
and charged lepton masses and mixings based on non-supersymmetric SO (10) without flavor 
symmetry was done The authors did not include the neutrino sector in the numirical 
fitting. Their result for atmospheric angle was sinOatm = 0.89. However, as this paper shows. 
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when the neutrino sector is included, the result is not only a better fit for atmospheric angle 
sin 6 atm = 0.776, but also the known light neutrino mass differences are accommodated. 

The paper is organized as follows. In section 2, a general structure of the fermion mass 
matrices for the second and third generations is constructed. Then, based on that structure, 
the fermion mass hierarchy and relations are explained. In section 3, it is shown that intro- 
ducing several 10-plets of matter fields to the model leads to the doubly lopsided structure 
which produces large neutrino mixing angles and small quark mixing angles simultaneously 
pT] . Then, some analytical expressions for quark masses and mixing angles at GUT scale are 
derived in a certain approximation to the model parameters. In section 4, exact numerical 
analysis is done to find the outputs at GUT scale. To get predictions of fermion masses 
and mixings at low scale, the quark masses and mixing at GUT scale will be run to the 
low scale by using renomalization equations. In section 5, It is shown how to get a suitable 
right-handed neutrino mass structure that gives the correct fitting for atmospheric angle 
after combining the charged lepton contribution. Section 6 is the conclusion. 

2 Fermion Mass Structure in SO(10) x A4 Symmetry 

In this section, the renormalizable Yukawa couplings of the SM fermions with the extra 
spinor-antispinor matter fields are considered as a concrete example of the model. The known 
matter fields of the SM (quarks and leptons) plus the right handed neutrino are contained 
in the three spinors (16j,3), where i stands for flavor index. The ordinary fermions, 16j, 
do not couple with 45^ in the minimal SO(IO). As a result, some of the predictions of 
minimal SO(IO) such as = rris and mc/rrit = rris/mb will follow, which are badly broken 
in nature. Therefore, extra heavy fermion fields must be introduced in order to allow the 45^^ 
to couple directly with the quarks and leptons of the standard model. The transformation 
of the ordinary fermions and the extra matter fields under A4 and the additional symmetry 
Z2 X Z4 X Z2 are summarized in Table 1. Let us consider first the invariant superpotential 
Wi under the assigned symmetry that contains the coupling of ordinary fermions with the 
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spinor-antispinor matter fields. 

Wi = 6il6iT6ilffi + 62l6iT62l'^i + J116iT6345jj + al63l62l0jf 

+M1I61T61 + M2I62T62 + M3I63T63. (1) 

Table 2 summarizes the transformation of the Higgs fields needed to achieve a minimum 
breaking scheme and the Higgs singlets needed to break symmetry. Although in this 
model, the structure in Eq.l does not include the Yukawa term 16il6il0if which is forbidden 
by the discrete symmetry Z2 x Z4 x Z2, the ordinary standard model fermions get their masses 
through their coupling with heavy extra fields. This is similar to how the light neutrinos 
get their masses through the coupling with the heavy right handed neutrinos in the known 
see-saw mechanism. 

The coupling terms in the superpotential Wi can be represented diagrammatically as 
shown in Fig.l. After integrating out the heavy states, the approximate effective operators 
can be read from the diagram. 

The VEVs of Higgs fields can be written down in general form as: 

{45h) = nQ, (3) 
ei 

(1^^) = I ^2 I , (4) 

(Ik) = I ^2 h (5) 

S3 

(5(10)) = Vu, (5(10)) = Vd. (6) 

Here the notation {p{q)) refers to a p of SU(5) contained in a g of SO(IO). The Q from Eq.3 
is a hnear combination of SO(IO) generators. One can redefine, without loss of generality, 
the hght fermion states as: 

16161 + 16262 + I6363 = 6I63 

I61S1 + I62S2 + I63S3 = S{16'2Se + IG'r^ce), (7) 
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Table 1: The transformation of the matter fields under SO(10)x744 and Z2 x Z4 x Z2. 



where e = y'ef + e| + e| and 5" = •\/sf+~sf+~s|. In terms of the redefined light fermion 
states, after dropping the prime notation and plugging in the VEVs, one gets 



Wn 



(163l62(5(163)Se + 163163(5(163)0^). 



(8) 



M1M2M3 

In general, the above effective operator can be written in terms of quark and lepton fields as 



Wp 



-{F^F^QpSe + F^F2QFcSe + F^F^{Qf + QFc)ce). 



(9) 



M1M2M3 

Here F is a general notation for up quarks {U), neutrinos (N), charged leptons (L) and down 
quarks {D). The quantity Qf{Qf<^) refers to the assigned charge of the left handed fermion 
(charge conjugate of the left handed fermions) after breaking SO(IO) group down to the the 
SM group. The unbroken charge Q can be written as a linear combination of two generators 
that commute with SUii)^ x SU{2)l x U{1)y as follows: 



6 Y 



(10) 



where I3R is the third generator of SU{2)ji and Y is the hypercharge of the Abehan U(l) 
group. The charge Q for different quarks and leptons is given below. 

2 



1 4 
Qu = Qd = -^S, Qu'^ = -1 - -5, Q^c = 1 + -5, 

3 6 
Qi^Qn^ -^(^, Qic ^1 + -S, Q„c = -1. 



Eq.lO can be expressed in the following matrix form: 
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Table 2: The transformation of the Higgs fields under SO(10)x744 and Z2 x Z4 x Zi. 



It is assumed that the masses of the extra heavy fermions are a little above the GUT scale, so 
the factors Np = that come from doing the algebra exactly are approximately 

equal to one, where T — VteS / MiM2M^. The first feature of the general mass matrix of 
the light fermions in Eq.l3 is an explanation for the mass hierarchy between the second and 
third generation in the limit sq ^ It is remarkable that a relation among generations is 
related to the vacuum alignment of A^ Higgs. 

Another feature of the above light fermion mass matrix m° = m° is obtained through 
M£)33 = Mi33, which follows from the relation Q^c -\- Qd — Qi'= + Qi- This relation occurs 
because both down quarks and charged leptons get their masses from the same Higgs. 

A further consequence of the light fermion mass structure is that m° 7^ m°. This inequal- 
ity relation follows from m°/m° = L32L23/ D22D23 = Qi'^Qi/Qd'^Qd which is not necessarily 
equal to one. This leads to the following question: What VEV direction should be given to 
45ij in order to obtain the Georgi-Jarlskog relation m° = 3m° ? There are two choices, either 
5 — > or 5 ^ —1.25. The former choice gives the unwanted relation (m°/m°)/(mg/m°) 1, 
while the later leads to {m^^/m^) / {m^g/m^) — >• 0. Thus, a good fit for 5 should be around 
-1.25. 

3 Extension to the First Generation and Doubly Lop- 
sided Structure 

In this section, the vector 10-plets fermions are added to the model to generate masses 
and mixings of the first generation. Adding these vector multiplets does not contribute to 
the up quark mass matrix since 10-plets do not contain a charge of (±2/3). Therefore, up 
quark matrix is still rank 2 and this is consistent with ^ pa 10~^ being much smaller than 
^ fa 10~^ and fa 0.3 x 10~^. First, it will be shown how the model leads to the doubly 
lopsided structure by employing these vector multiplets, then some analytical expressions 
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for masses and mixing angles of fermions at GUT scale will be derived. Let us first consider 
the invariant couplings under the assigned symmetry which can be read from the Feynman 
diagram in Fig. 2. The allowed couphngs in the superpotential W2 are 

W2 = leaoae^f + MiolOaO^ + hfijk^O'iiO'jlHk + hijklOilOjlHk. (13) 

The important point is that Fig. 2 gives a fiavor-symmetric contribution to the down quark 
and charged lepton mass matrices. In order to understand this, recall that the general 
product of three triplets: (oi, 02, 03), (61, &2, ^'3), and (ci, C2, C3) that transform as a singlet 
under A4 is given by 

hi{a2bsci + 036102 + 016203) + h2{a3b2Ci + 0163C2 + 0261C3). (14) 

The third term of Eq.l3 gives a symmetric contribution since there are two identical 10- 
plets. The last term in Eq.l3 has been ignored by assuming the Yukawa couplings hijk are 
very small. The contribution of Fig. 2 to the mass matrices of the down quarks and charged 
leptons after integrating out the extra vector multiplets is then 

/ C12 Ci3 \ 

M£ = Ml, oc C12 C23 , (15) 

\ Ci3 C23 / 

where C12, C13, and C23 are proportional to ei, 62, 63 respectively. To obtain the desired 
fermion mass structure (the doubly lopsided structure, which is going to be explained later 
in this section), other couplings are needed to be included by employing the four vector 
10-plets plus adding another Higgs singlet I'^jj to the model (their transformations under 
the assigned symmetry are shown in Tables 1 and 2). The purpose of these couplings is 
to give a flavor-antisymmetric contribution to the down quarks and charged leptons mass 
matrices. Since the adjoint of SO(IO) (45//) is an antisymmetric tensor which changes the 
sign under the interchange 10- ^ 10"'. One can consider employing the Yukawa coupling 
lOf 10^45//. Also, from the fact that when we write the SO(10)-vectors in the SU(5) basis 
such as lOj = 5i + 5i, the charged lepton and down quark content of 5^ or 5j have different 
chirahty, therefore the structures of matrices Ml and Mo have opposite signs (look at the 
mass structures in Eqs. 18-19). It is important to emphasize that the minimum Higgs breaking 
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scheme assumption does not allow us to add another adjoint to the model. Therefore, the 
same adjoint 45 h Higgs representation that breaks SO (10) group to the SM group is going 
to be used. Additional couplings to the previous superpotential can be read from Fig.3. 

Ws = io;io;'i'^, + mio^'ior + 10^10:45^, (le) 

where (45ij) has been defined previously. The following VEV to the Higgs singlet 1'^ is given 
below: 

m = M2 I . (17) 
S3 

After integrating out the heavy states, the following contribution to the Ml and is 
obtained: 

-6sQi S2Q1 \ 
Mt oc ( 5^Qi -5^Qi , (18) 
-52Q1 SiQi / 

oc I -SsQd^ 6iQd^ , (19) 
S2Qd<= -SiQd'= / 

where the overall constant has been absorbed in the redefinition of 5i, 62 and ^3. Eqs.18- 
19 show that the off-diagonal elements of (Mf) are proportional to Q^c (Qi)- This is 
because 5j(10) contains in its representation the charge conjugation of a color triplet of the 
left handed down quarks d% and the left handed charged leptons eu- The full tree level 
mass matrices, which are obtained by adding the three superpotentials Wi + W2 + W3, have 
the following forms: 

Cu + SS^i^) -52« + C 

Mi = m°| Ci2-353(^) 5ia-3s{^)+P | , (20) 

(-62^ si^)+6i{^)+P 1 

c,2 + Ss{^) -2<52(^) + C 

Mn ^ m",{ cu-Ss{^) 25,{^) + s{^) + /3 |, (21) 

C 1 



Mu ^ mil C-^)s I , (22) 

{^)s 1 



Mjv = m° I (^^)s I , (23) 
s 1 



where the convention being used here is the left handed fermions multiphed from the right. 
The parameters of the model have been defined as follows: 

5 = -l + a, (24) 
(|5 + l)c,- 

The above fermion mass structures have eight parameters. If a goes to zero, the fermion 
mass matrices in Eqs. 20-21 go to the SU(5) limit {m^ = m°, = m°, m[J = m'^). To 
avoid the bad prediction of SU(5) for lighter generations, the good numerical fitting for a 
should deviate from zero. On the other hand, to keep the good SU(5) prediction for third 
generation, the parameter a should satisfy a « 1. li 5i and 82 are of order one and the other 
model parameters are very small C, a, ^3, cu, s « Si, S2), the model leads to the doubly 
lopsided structure. To see this clearly, let us go to the limit where the small parameters are 
zero (except s) . So the Mjj and Ml go to the following form: 

/ \ 

Mi = = m° (f ) . (25) 

(f)+5i(|) 1 J 

In diagonahzing Ml of the Eq.25, the large off diagonal elements Si and S2 that appear 
asymmetrically in M^ and Ml must be eliminated from the right by a large left handed 
rotation angle Ogoi in the 1-2 plane, where ta'R{9soi) = The next step of diagonalization 
is to remove the large element o" ~ (5f + 5|)^ that has been produced after doing the first 
diagonalization where (3,2) element of the matrix in Eq.25 is replaced by a. This can be done 
by a rotation acting from the right by a large left handed angle $23 in the 2-3 plane, where 
tan(^23) ~ — c- On the other hand, there are no corresponding large left handed rotation 
angles in diagonahzing Md since Ml = M^. However, the large off diagonal elements in 
Ml) can be eliminated by large right handed rotation angles acting from the left on the M^ 
in Eq.25, while the left handed rotation angles are small. This explains how the doubly 
lopsided structure leads to the small CKM mixing angles and the large neutrino mixing 
angles simultaneously. If the parameters C12, ^s, and ( are zero, analytical expressions can 
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be written down for the ratios of quark and lepton masses of the second and third generations, 
Vcb, and neutrino mixing angles ( tan and tan ^23) in terms of Si, S2, s ,a, and /3: 

_ g^(l-Q;)(l + 4a) 
~ 25 ' 



ml 5(l + A(3 + 2aW + 5i)) ' 



/5 + 



25 

s(3+2a) 



:i + A(3+2«)2(5f+5i))' 

-s(l + 4q;) 



tan^ 



52( 



6— a ^ 



6 — a. 



12 = 

tan^23 = + 

These expressions are derived by using the approximation a,s,/3 << 81,82 and are useful 
for fitting the data. The best fit for the data is obtained by setting tan ^23 = —2 and 
tan ^12 = 0.68, which correspond to ^23 = —63° and Qyi — 34". The central value of the 
atmospheric angle is around 45°. In order to bring 63° down close to the central value, the 
neutrino sector is required to be included as shown in section 5. Also, it will be shown that 
the contribution of the neutrino sector to the solar angle is small. 

4 The Numerical Results 

The model can be shown to be concrete by giving numerical values to the parameters of the 
model, and producing the six mass ratios of quarks and leptons, CKM mixing angles (V^s, 
Vuhi and V^b), CP violation parameter 77 = —IineVubVcs/VusVcb), and neutrino mixing angles 
(sin ^12, and sin ^13). The fermion mass texture in Eqs. 20-23 has eight parameters {5i, 62, Ss, 
a, j3, s, C, and C12). These parameters are going to be used to fit 12 quantities. If 61 = —1.302, 
S2 = 1.0142, 53 = 0.015 X e^-95\ a = -0.05801, s = 0.29, ( = 0.0105, C12 = -0.00153e^-^^26j^ 
and /3 = -0.12303, the following excellent fit at GUT scale is obtained : ^ = 0.002717, 
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Table 3: This Table shows the comparison of the model predictions at low scale and the 
experimental data. 

4 = 0.958, 4 = 0.00473, ^ = 0.0585, ^ = 3.63, 4 = 0.302, r] = 0.357, = 0.2264, 
Kb = 0.0037, Vch = 0.0362, sin^ia = 0.569, and sin^ig = 0.0653. The above numerical 
fittings lead to sin = 0.904 which is not close to the central value sin^^23™ = 0.707. 
One can see from the superscript L that the mixing angle ^23 comes only from the charged 
lepton contribution. To obtain close to the expected atmospheric angle and the correct 
neutrino mass differences, it is important to include the neutrino sector contribution to the 
atmospheric angle by finding out a suitable right handed neutrino structure which respects 
the assigned symmetry of the model. 

In order to compare with experiment, the predicted fermion masses and mixing angles 
at the low energy scale are needed to be found. The above numerical values of the fermion 
masses and mixing angles which are obtained at GUT scale have been evolved to the low 
scale in two steps. First, the running from GUT scale to Msusy — 1 TeV is done by using the 
2-loop MSSM beta function. The running factors denoted by rji depend on the value of tan j3. 
The known fermion masses and mixing data are best fitted with tan[/3] = 10. The running 
factors for tan[/3] = 10 are (773/6, Vh/t, Vb/r: Vc/t: Vcb— rjub)—{0.S736, 0.9968, 0.5207, 0.73986, 
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0.910335), where r]i/j = (m°/m°)/(mj(iret;)/mj(lTet;)) and r]cb,ub = V^f^,^jVcb,ub{iTev). 
The second step is to evolve the fermion masses and mixing angles from Msusy = 1 TeV to 
the low scale. The renormalization factors rji that run fermion masses from their respective 
masses up to the supersymmetric scale Msusy = 1 TeV are computed using 3-loop QCD and 
1-loop QED, or electroweak renormalization group equation (RGE) with inputs as{Mz) = 
0.118, a{Mz) = 1/127.9 and sm6w{Mz) = 0.2315. The relevant renormalization equations 
can be found in [12] [IB]- The results are {rjc, rjb, rj^, rj^, 77^, rjt, ■?7^,6=?7cb) = (0.4456, 0.5309, 
0.8188, 0.83606, 0.8454, 0.98833, 1.0151). 

By using the above renormalization factors, = 1776 MeV, and rrit = 172.5 GeV, the 
following predictions at low scale can be obtained: mc{m^ = 1.4 GeV, mb{mb) = 5.2 GeV, 
me(me) = 0.511 MeV, m^(m^) = 105.6 MeV, md{2 GeV) = 7.5 MeV, m,(2 GeV) = 132 MeV, 
ri = 0.357, Vus = 0.2264, Kb = 0.004, Vcb = 0.0392, sin^ia = 0.569, and sin^ia = 0.0653. 

Note that the numerical value of rrib is not in agreement with the experimental value 
rUb = 4.20lo'j7 GeV [14J. In order to fix this, the finite gluino and chargino loop corrections 
[15] are required to be included to the down type quark masses ( ma, rris, riib ), which are 
denoted respectively by {l+Aa), (1+As), (l+A;,). These corrections are proportional to 
the supersymmetric particle spectrum: A5 ^ tan/9(^^^2^^^p— [/(m? /M|) — /(m? /M?)] + 

16^ m^^-m^ [/(m^ / fi^) — f {m'^ //^^)])) where /(x) = ln(x)/(l— a;) and the first (second) term 
refers to gluino (chargino) correction. Similar expressions exist for and A^, but without 
the chargino contribution and b ^ s,d. If the chargino loop corrections are negligible and mj, 
nis, and are degenerate, the equality relation A^ = A^ = A;, is approximately satisfied. 
In order to get better fitting for down type quark masses, let us take Ad = As = Ab = —0.17 
which gives 771^(2 GeV) = 6.24 MeV, m°(2 GeV) = 109.65 MeV, and mb{mb) = 4.31 GeV. 
The comparison of the model predictions and experimental data at low scale is summarized in 
Table 3, where the quark and charged lepton masses, the CKM mixing angles (Vub, Ks, Vcb), 
the neutrino mixing angles [smOgoi, sinOatm, sin 6*13), and the CP violation parameter (77) are 
taken from [T3]. The masses are all in GeV. Although the model here predicts my_{GUT) = 0, 
the quantity mud = ("^u+"^d)/2 is considered in Table 3, where it is assumed that the tiny up 
quark mass at GUT scale may be generated either by including the coupling 16il6jl0// into 
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the model or by considering higher dimensional operators. Suppose m„(2 GeV) = 2.4 MeV, 
the model predictions of the quantities rriud and which are well-known from lattice 
calculation [12], are given in Table 3. The asterisks in Table 3 indicate that the model 
predictions of neutrino mixing angles are obtained after including the neutrino sector in 
section 5. 

5 Right Handed Neutrino Mass Structure 

Up to now, the model gives excellent agreement with the known values for the CKM mixings, 
the quark masses, the charged lepton masses, CP violation parameter, and the neutrino 
mixing angles (sin5i2 and sin ^13). However, the whole picture is still not complete and 
the following question arises. What is the appropriate light neutrino mass matrix [M^, = 
-MJ^M^^Mn) that gives not only the correct contribution to the atmospheric angle, but 
also the correct neutrino mass differences: Am^2 = (7.59 ±0.2) x 10~^ eV"^, ^^^32 = (2.42 ± 
0.13) X 10~^ eV"^ [11]? In the other words, we are looking for a suitable structure of right 
handed neutrino mass matrix Mr since is fixed. Recall that MNS mixing matrix is given 
by 

Umns = UlU,, (27) 

where Ul and Ui, are the unitary matrices needed to diagonalize the hermitian lepton matrix 
mI^Ml and the light neutrino matrix M^^ respectively. 

^diag^^diag ^ uImIMlUl, Mf""^ = U^M^u, (28) 

where Mj, is assumed to be real and symmetric. The Dirac neutrino mass matrix Mj^ in 
Eq.23 has vanishing first row and column and the same is true for M^,. So the matrix required 
to diagonalize M^, is simply a rotation in the 2-3 plane by an angle 6^, while Uj^ is determined 
numerically from the charged lepton mass matrix. Thus, the mixing matrix of neutrinos is 
given by 

/ -0.1437 -0.80681 0.1353 + 0.55301 0.0653 \ / I \ 

Umns=\ 0.2497 + 0.05821 0.3394 -0.04001 0.9041 cos 61^ sin 6*^ . (29) 

V -0.5125 + 0.00021 -0.7477 + 0.000061 0.4222 / \ -sin 6*^ cosOy ) 

One can conclude the correct contribution of the neutrino sector to the atmospheric angle 

is around 6y=—2Q°. For example, if we take 6*1^=— 20° the neutrino mixing angles (sin 6'atm, 
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sinOsoi, sm6'i2) become (0.707, 0.53, 0.21). In order to find the suitable right handed neutrino 
mass structure, one can easily prove the inverse of the see-saw relation. 

Mr = -MNU,{Mf-s)-'UjMj. (30) 

A similar technique was used in ref [IT]. Note that one of the eigenvalues of M^, is zero (i.e. 
]\^dzag jg singular), so the inverse of M^''"-^ does not exist. To overcome this problem, one 
can generally define M^*"^=diag( mi, m2, ) and mi will not appear in M/j. By using the 
numerical result of Mjv, 6'^=— 20°, and m2/m3=0.178, the right handed mass structure can 
be presented numerically. 

\ 

0.0186 -0.13 . (31) 
-0.13 1 / 

From the above numerical mass matrix, one concludes {Mji)23 x (M/j)23 ~ (^^2)22, so to 
a good approximation, the above numerical structure can be represented analytically as 
follows: 

\ 

r2 ar . (32) 
ar 1 / 

The constant a should not be equal to one because then Mr would be singular, but the 
constant a should be around 1. Now our mission is to find the Yukawa couplings that 
respect the sjTiimetry of the model and lead to an analytical structure similar to Eq.32. 
This can be accomplished by considering the following Yukawa couplings represented by 
Feynmann diagram in Fig.4. 

W4 = 16il6Hli + hijkhl-lHk + (33) 

Where two fermion singlets Ij and 1^ which couple with the singlet Higgs 1-'^ have been 
introduced (their transformation under SO(10)xy44 and the additional symmetry are shown 
in Tables 1 and 2). The singlet transformation of the product of three triplet under A4 of 
the second term of Eq. 33 is given by hi^NiN^a^ + N2N^ai + N-iN[a2) + /i2(A^iA^|a2 + 
A^sAi'lai + N2Nla3), where ai, 0^2, and as are the VEV's components of 1"^. By assuming 
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/ii=/i2, Fig.4 leads to the desired right handed neutrino mass structure. 



Mr^K 



( "1 n-,n-..r-l I .. ^. .. 'I -«iK-aj+«i) \ 

/-I I 2 N a| -02(-Qi+ai+a3) 

I -Qi("i-Qf+"3) -a2(-af+"2+"3) ]_ J 

\ «3(Q^+a^+a§) a3(a^+a|+a|) / 



(34) 



By comparing the 2-3 block of the above structure with the mass structure in Eq. 32, one 
can see the constant a is equivalent to the quantity ((— + q;| + ctg) / (a^ + q;| + ctg)) which 
is equal to one in the limit a\ 0. So, let us expand the eigenvalues of the right handed 
neutrino mass structure in Eq.34 around a\. 



a 



2 ' ^,2/' ^,2 I ^,2^2 ' ^V-^iy; 

3 



- RT^+(ai + air/2 

One can see two of the right handed neutrino masses are approximately degenerate for small 
values of CKi (i.e. Mr2 ^ Mrs). By setting (ai, q;2, ctg, A)=(-0.05, 0.125, 0.994, 8.42 x 10^^), 
the numerical fit for the neutrino mixing angles, the light neutrino masses and right handed 
neutrino masses are obtained as follows: 



mi 


= OeV, 


sin 6soi 


= 0.551, 


Mri 


= 8.57 X 10^^ GeV, 


1712 


= 0.01 eV, 


sin 9atm 


= 0.776, 


Mr2 


= 1.3 X 10^2 GeV, 


1713 


= 0.056 eV, 


sin 


= 0.154, 


Mr3 


= 1.28 X 10^^ GeV. 



(36) 



6 Conclusion 



Both largeness of neutrino mixing angles and the smallness of CKM mixing angles can be 
accounted for by a lopsided structure of the charged leptons and down quarks. This structure 
was obtained in several studies and has been used for many models of fermion masses. The 
model that has been studied here is the first SUSY grand unification model based on the 
gauge symmetry SO (10) with the discrete family symmetry A4 that leads to the doubly 
lopsided structure. A few works on SO(IO) XA4 have recently been pubhshed, but what 
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makes this work unique is the assumption of using the minimal set of Higgs fields that 
break SO(IO) to the SM group. This assumption acts as an important guide for searching 
for good models. The possibilities of renormalizable Yukawa interactions for quarks and 
leptons are very limited because the minimum Higgs breaking scheme is imposed and the 
superpotential must respect the assigned symmetry of the model. Based on that, a general 
mass structure for the heavy generations has been obtained which explains the following 
features: (1) m° ~ m°, (2) ^ = 3, (3) ^ << It is important to mention that another 
work [2] got the same mass structure for heavy fermions. In that work, the authors did 
not employ the flavor symmetry and showed that the hierarchy between the second and 
third generations can be understood by choosing the specific direction of (45/^). Also, they 
employed another adjoint Higgs field (45//) to include the first family to their model. On the 
other hand, in this presented study, the above features of heavy fermions have been obtained 
by picking specific direction of (45//), but the hierarchy between the three generations can 
be understood in the framework of y44-symmetry. Without adding another adjoint to the 
model, the first family is successfully included to the model and excellent predictions are 
obtained. 

For fitting purposes, some approximate analytical expressions given in Eq. 26 are derived 
for mass ratios and mixing angles of the quarks and the leptons by combining the Yukawa 
couplings represented by three Feynman diagrams in Figures 1-3. However, exact numerical 
fitting at low scale was done. Without including the neutrino sector, the model predictions at 
low scale for the masses and the mixing angles (except the atmospheric angle) of the quarks 
and the charged leptons, as well as CP violation parameter are in excellent agreement with 
data (i.e. within 2cr). The atmospheric angle needes to be corrected by considering the 
neutrino sector. The symmetry of the model succeeds to produce the appropriate right 
handed neutrino structure that gives not only the correct contribution to the atmospheric 
angle, but also the correct neutrino mass differences. The neutrino contribution to the solar 
angel is negligible. Besides, after combining the two contributions (neutrinos and charged 
leptons contributions) the neutrino mixing angle 9is is predicted. 
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Figure 1: This figure shows diagrammatically the couphngs in the superpotential Wi 
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Figure 2: This figure leads to the fiavor symmetric contribution to the down quarks and 
charged leptons. 
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Figure 3: This figure leads to the fiavor antisymmetric contribution to the down quarks and 
charged leptons. 
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Figure 4: This figure leads to the right handed neutrino mass matrix. 
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